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Increased serum uric acid (SUA) is a risk factor for gout and renal and cardiovascular 
disease (CVD). The purpose of this study was to identify genetic factors that affect 
the variation in SUA in 632 Mexican Americans participants of the San Antonio Family 
Heart Study (SAFHS). A genome-wide association (GWA) analysis was performed using 
the lllumina Human Hap 550K single nucleotide polymorphism (SNP) microarray. We 
used a linear regression-based association test under an additive model of allelic 
effect, while accounting for non-independence among family members via a kinship 
variance component. All analyses were performed in the software package SOLAR. 
SNPs rs6832439, rs13131257, and rs737267 in solute carrier protein 2 family, member 
9 {SLC2A9) were associated with SUA at genome-wide significance (p< 1.3 x 10^''). 
The minor alleles of these SNPs had frequencies of 36.2, 36.2, and 38.2%, respectively, 
and were associated with decreasing SUA levels. All of these SNPs were located in 
introns 3-7 of SLC2A9, the location of the previously reported associations in European 
populations. When analyzed for association with cardiovascular-renal disease risk factors, 
conditional on SLC2A9 SNPs strongly associated with SUA, significant associations were 
found for SLC2A9 SNPs with BMI, body weight, and waist circumference (p < 1 .4 x lO^'^) 
and suggestive associations with albumin-creatinine ratio and total antioxidant status 
(TAS). The SLC2A9 gene encodes an urate transporter that has considerable influence 
on variation in SUA. In addition to the primary association locus, suggestive evidence 
(p < 1 .9 X 10^^) for joint linkage/association (JLA) was found at a previously-reported urate 
quantitative trait locus (Logarithm of odds score = 3.6) on 3p26.3. In summary, our GWAS 
extends and confirms the association of SLC2A9 with SUA for the first time in a Mexican 
American cohort and also shows for the first time its association with cardiovascular-renal 
disease risk factors. 



Keywords: variance components decomposition approach, joint linl<age/association analysis, l<inship, 
hyperuricemia 



INTRODUCTION 

Hyperuricemia is a risk factor for gout, renal disease, and cardio- 
vascular disease (CVD) (Cirillo et al., 2006; Nakagawa et al., 2006) 
and is known to aggregate in families (Dbcon, 1960; Friedlander 
et al., 1988; Cameron and Simmonds, 2005). Variation in serum 
uric acid (SUA) levels is controlled by both genetic and envi- 
ronmental factors. Family-based studies have reported significant 
heritabilities for SUA levels with estimates ranging from 25 to 
73% (Rao et al, 1982; Rice et al, 1990; Wilk et al, 2000; Tang 
et al, 2003; Yang et al, 2005; Nath et al., 2007; Voruganti et al., 
2009a,b). 

Genome-wide studies conducted to identif}' significant 
linkages for the variation in SUA have found quantitative trait 



loci (QTL) on several chromosomes for different populations; 
a QTL on chromosome 15 in Framingham Heart Study (Yang 
et al., 2005); on chromosome 8 in the Genetic Epidemiology 
Network of Arteriopathy (GENOA) (Rule et al., 2009); on chro- 
mosome 4pl5 in an Australian Cohort (Cummings et al., 2010). 
However, none of them were in the same location as found 
previously in this Mexican American Cohort (Voruganti et al., 
2009a). Similarly, genome-wide association (GWA) and candi- 
date gene studies have found several genes to be associated 
with SUA, mainly solute carrier protein 2 family, member 9 
(SLC2A9) (Li et al, 2007; Brandstatter et al., 2008; Caulfield 
et al, 2008; Dehghan et al, 2008; Doring et al, 2008; McArdle 
et al., 2008; Stark et al, 2008; Vitart et al, 2008; Wallace et al. 
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2008; Kolz et al, 2009; Charles et al, 2011; Karns et al., 2012), 
ATP-binding cassette, subfamily G, member 2 {ABCG2) (Kolz 
et al, 2009; Woodward et al, 2009; Nakayama et al, 2011; 
Karns et al., 2012; Zhang et al., 2013); solute carrier protein 
22 family, members 11 and 12 {SLC22A11, SLC22A12) (Kolz 
et al., 2009; Tin et al., 2011), solute carrier protein 16 family, 
member 9 (SLC16A9), PDZ domain containing 1{PDZK1) (Kolz 
et al., 2009). A large GWAS in another European cohort, Kottgen 
et al. (2013), confirmed or replicated 28 previously reported 
associations for SUA. However, most of these studies were con- 
ducted in Caucasian, African American, or Asian populations. To 
date no GWAS of SUA levels have been conducted in Mexican 
Americans — a population affected with disproportionate bur- 
den of CVD and type 2 diabetes (Flegal et al, 1991; Go, 2013; 
http://www.diabetes.org/diabetes-basics/diabetes-statistics/). 

SUA is a strong marker of risk for renal disease and is consid- 
ered clinically relevant and of high significance. Its role as a major 
risk factor in the development of kidney stones or nephrolithiasis 
and gout is well-recognized. The relationship between uric acid 
and kidney function seems to be two-sided. On one hand, decline 
in glomerular filtration rate (GFR) (kidney function parameter) 
may lead to elevation of uric acid, on the other hand, increase in 
uric acid seems to alter glomerular function through renal vaso- 
constriction and increased rennin expression (Sanchez-Lozada 
et al., 2002). Hyperuricemia is also known to induce endothelial 
dysfunction and inflammation indicating a role in atherosclero- 
sis as well as chronic kidney disease (Ejaz et al., 2012). Rodent 
Sanchez-Lozada et al, 2002; Ejaz et al., 2012 and human stud- 
ies (Nakagawa et al., 2005; Liebman et al., 2007) have shown the 
role of hyperuricemia in hypertension, atherosclerosis, CVD, ini- 
tiation and progression of renal disease, and metabolic syndrome. 
In fact, drugs that decrease SUA levels, mainly aUopurinol, have 
been shown to improve survival in chronic heart failure patients 
(Gotsman et al., 2012), improve endothelial function in patients 
with chronic kidney disease (Yelken et al, 2012), and reduce 
oxidative stress and improve endothelial function in patients with 
coronary artery disease (Rajendra et al., 2011). Elaborating on the 
protective effects of reduced SUA levels on end-stage renal disease 
(ESRD), Doria and Krolewski (2011) point out that lowering of 
SUA may be a potential approach to treat ESRD in patients with 
diabetes. 

Here we present results of the first GWAS in Mexican 
Americans of the San Antonio Family Heart Study (SAFHS) that 
attempts to identify genetic factors that affect the variation in 
SUA. We extend these findings to the analysis of association of 
SLC2A9 with cardiovascular and renal factors given the role of 
SLC2A9 in hypertension and renal urate transport. 

MATERIALS AND METHODS 

Population characteristics: The San Antonio Family Heart Study 
(SAFHS) was initiated in 1991 to identify genes influencing 
the risk of CVD in Mexican Americans. Study subjects have 
been recalled up to three times to acquire longitudinal data; 
the analysis reported here is based on the second recall, 2002- 
2006. Individuals were recruited from large Mexican American 
families, residing in San Antonio, TX without regard to dis- 
ease status. Participants in this study were recruited from 40 



extended families, with probands between 40 and 60 years of age 
(MacCluer et al., 1999; Mitchell et al, 1999). Eligibility criteria 
required the proband to have at least six first-degree relatives 
(excluding their parents) 16 years or older and who resided 
in San Antonio, TX. At each recruitment phase, subjects were 
brought to a research clinic at the University of Texas Health 
Science Center — San Antonio (UTHSCSA) for interview and 
examination by trained recruiters and nurses. Anthropometrics 
including height, weight, and waist circumference; blood pres- 
sure; and self-reported information regarding medical history 
and socio-demographic status were obtained in all phases of 
SAFHS data collection. Blood samples were collected from all 
participants after an overnight fast and plasma and serum were 
prepared and stored at — 80°C until analyzed. Blood samples 
were also drawn at 2 h after a standard oral glucose toler- 
ance test; diabetes was diagnosed if the 2h glucose level was 
ll.lmmol/1 or higher, or if the subject had been prescribed 
antidiabetic medication. The final analysis sample with complete 
phenotype and genotype data included 632 SAFHS participants. 
Written informed consent was obtained from all participants 
to participate in this study. AH research and consenting pro- 
tocols were approved by the Institutional Review Board of the 
UTHSCSA. 

PHENOTYPING 

Uric acid was measured in serum by a colorimetric assay using 
uricase and peroxidase (Domagk and Schlicke, 1968). Serum 
creatinine was estimated by the modified kinetic Jaffe reaction 
(Beckman Synchron LX System). GFR was estimated by the 
MDRD equation: eCFR (ml/min/1.73 m2 body surface area) = 
186 X serum creatinine x age x (—0.203) x (0.742 if female) x 
(1.210 if black) (Arar et al., 2008). Cardiovascular risk factors 
included blood pressure, weight, waist circumference as well as 
fasting plasma levels of glucose, insulin, lipids, and cytokines 
measured by standardized reference procedures (Arar et al, 2008; 
Voruganti et al, 2009a). A single-void morning urine sample 
was collected from each participant for measuring albumin and 
creatinine and estimating albumin to creatinine ratio (UACR). 
Description of their measurement techniques are given in Arar 
et al. (2008). Indicator variables were coded for diabetes diagno- 
sis and for data from self-report and medical history on current 
use of blood pressure medication, alcohol consumption, and 
smoking. Covariates of SUA were included in analysis models as 
described below and in the Results. 

GENOME-WIDE ASSOCIATION (GWA) ANALYSIS 

GWA analysis was conducted in the SAFHS based on SNP geno- 
types obtained using the lUumina HumanHapSSO BeadChip 
(lUumina, San Diego, CA). Our experimental error rate (based on 
duplicates) was 2 per 100,000 genotypes. The average call rate per 
individual sample was 97%. Approximately 1 per 1000 genotypes 
was blanked due to Mendelian errors. Specific SNPs were removed 
from analysis if they had call rates <95% (about 4000 SNPs) or 
deviated from Hardy-Weinberg equilibrium at 5% false discov- 
ery rate (FDR) (12 SNPs). Missing genotypes were imputed from 
pedigree data using MERLIN (Abecasis et al, 2002). SNP geno- 
types were checked for Mendelian consistency using the program 
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SimWalk2(Sobel and Lange, 1996). The estimates of the allele fre- 
quencies and their standard errors were obtained using SOLAR 
(Almasy and Blangero, 1998). 

MEASURED GENOTYPE ANALYSIS 

Each SNP genotype was converted in MERLIN (Abecasis et al., 
2002) to a covariate measure equal to 0, 1, or 2 copies of the 
minor allele (or, for missing genotypes, the weighted covari- 
ate based on imputation). These covariates were included in the 
variance-components mixed models for measured genotype anal- 
yses (MGA; Boerwinkle et al., 1986) vs. null models that incorpo- 
rated the random effect of kinship. For the initial GWA screen, 
SUA levels were regressed on selected covariates (see Results) 
and the association of the inverse-normalized residuals with each 
SNP covariate was tested independently as a 1 degree of freedom 
likelihood ratio test. Empirical thresholds for genome-wide sig- 
nificant and suggestive evidence of association were based on the 
distribution of p-values from 10,000 simulated null GWAS (i.e., 
simulations of a heritable trait with no modeled SNP covariate 
effects using the SAFHS pedigree and genotypes). The threshold 
for significance (p < 1.3 x 10^'') was defined as the cutoff for the 
lower 5% tails of the empirical distribution, and the threshold for 
suggestive evidence (p < 1.6 x 10^^) was the minimum p-value 
obtained not more than once per genome scan. 

LINKAGE AND JOINT LINKAGE/ASSOCIATION ANALYSIS 

Multipoint linkage analysis (Almasy and Blangero, 1998) was 
performed in SOLAR using estimates of locus-specific allele shar- 
ing based on genotypes for 461 STR markers. In addition, we 
employed a novel joint linkage/association (JLA) analysis for 
each SNP that tested each saturated model (including linkage 
and the fixed effect of the SNP) against a null model in which 
both effects were constrained to zero. Regression parameters for 
selected covariates (see Results) were estimated simultaneously 
with the linkage and association parameters. Because the linkage 
variance parameter was tested on its lower boundary, the dis- 
tribution of the likelihood ratio test statistic is distributed as a 
1:1 mixture of chi-square distributions with 1 and 2 degrees of 
freedom, respectively (Self and Liang, 1987). 

ANALYSIS OF CARDIOVASCULAR AND RENAL DISEASE RISK FACTORS 
FOR ASSOCIATION WITH SLC2A9 SNPs 

Cardiovascular and renal disease risk factors that were included in 
the analysis were anthropometric measures such as body weight, 
BMI, waist circumference, circulating levels of lipids such as 
triglycerides, high, low-density lipoprotein and total cholesterol, 
glucose, insulin, and kidney function phenotypes (serum crea- 
tinine, eGFR, and UACR). For these analyses, the appropriate 
significance level was determined to be 1.4 x 10^^ based on the 
number of SNPs investigated in SLC2A9. This significance value 
was computed taking into account the linkage disequilibrium 
(LD) pattern of these SNPs (Figure 1). 

RESULTS 

The mean ± SE age and SUA levels of participating individuals 
(m = 632) were 47.87 ± 14.8 yrs and 5.35 ± 1.38 mg/dl, respec- 
tively, with men having higher levels of SUA than women (6.1 ± 
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FIGURE 1 I Linkage disequilibrium pattern of SLC2A9 SNPs 
investigated in thiis study. 



1.7 vs. 4.94 ±1.6 mg/dl). Significant heritability was detected for 
SUA levels (h^ = 0.39, p = 2.3 x 10"*^) with age, sex, age^, inter- 
action between age and sex, waist circumference, systolic blood 
pressure, plasma triglyceride, plasma high-density lipoprotein, 
serum creatinine, and indicator variables for use of blood pres- 
sure medication, alcohol use, smoking status, and diabetic status 
as covariates. 

GENOME-WIDE LINKAGE/ASSOCIATION ANALYSIS 

Association and JLA tests were performed using a custom 
script for SOLAR (Almasy and Blangero, 1998) as described 
in methods. AH models included the random kinship as well 
as fixed effects of SNPs. Prior to analysis, SUA levels were 
regressed on age, sex, age^, interaction between age and sex, 
waist circumference, systolic blood pressure, plasma triglyc- 
eride, plasma high-density lipoprotein, serum creatinine, blood 
pressure medication, alcohol intake, smoking status and dia- 
betic status, and the residuals were inverse-normal trans- 
formed. The distribution of p-values from GWAS of SUA 
showed no evidence of inflation due to population stratification 
(Figure 2). 

SNPs rs6832439, rsl3131257, and rs737267 in SLC2A9 were 
associated with SUA levels at genome-wide significance (p < 
1.3 X 10"^) (Table 1, Figures 3, 4). The minor alleles of these 
SNPs had frequencies of 0.36, 0.36, and 0.38, respectively, and 
were associated with decreasing SUA levels. One additional SNP 
rs6449213 showed suggestive association (p < 1.6 x 10^*) with 
SUA levels. All these SNPs are located in introns 3-7 of SLC2A9. 
The allele frequencies of SNPs that had significant association 
with SUA are shown in Table 1. Minor alleles of the main 
four SNPs were associated with lower SUA levels (Table 2). We 
had previously reported significant evidence of linkage for SUA 
(LOD = 3.5) on chromosome 3p (Voruganti et al, 2009a). Using 
the JLA test, there was suggestive evidence of association of SUA 
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with SNPs in CNTN4, a gene within the 1-LOD support interval 
of our peak linkage signal (Table 1, Figure 5). 

ASSOCIATION OF SLC2A9 SNPs WITH CARDIOVASCULAR AND RENAL 
DISEASE RISK FACTORS 

In addition to SUA, we also tested the association of SLC2A9 
SNPs with several renal and CVD-related risk factors such as 
anthropometrics, glucose, insulin, lipids, cytokines, and renal 
function phenotypes. Of these, SNPs in the SLC2A9 gene were 
significantly associated with body weight, BMI, waist circum- 
ference, total antioxidant status (TAS) and urinary albumin- 
creatinine ratio (UACR) (Table 3). However, the SNPs associated 
with SUA levels were not the same as those associated with 
body weight, BMI, waist circumference, or ACR. Interestingly, 
the two SNPs associated with TAS were also associated with SUA 
levels. 




0 5 10 15 20 

Chi square - Quantiles 
(Inflation factor = 0.998) 



FIGURE 2 I A Q-Q plot showing the absence of inflation due to 
population stratification. 



When association was conducted, conditional on SNPs 
associated with SUA, SLC2A9 SNPs showed evidence of signifi- 
cant association with BMI, waist circumference and body weight, 
and suggestive association with TAS and UACR (data not shown). 
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FIGURE 3 I Genome-wide association analysis of serum uric acid. 
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FIGURE 4 I Association of serum uric acid with SLC2A9 SNPs. 



Table 1 | Joint linkage-association analysis of serum uric acid (Empirical genome-wide significance: p < 1.3 x 10 ^). 



Gene SNP° Coordinates (bp) MGA'' (p-value) JLA'^(p-value) Minor allele/frequency 



SLC2ASf' 


rs6832439 


9924319 


6.0 X 10-9 


2.7 X 10-8 


A/0.36 


SLC2A9 


rsl3131257 


9981889 


2.0 X lO"'' 


8.1 X 10-8 


T/0.36 


SLC2A9 


rs737267 


9934744 


4.2 X 10"^ 


1.7 X 10-7 


A/0.38 


SLC2A9 


rs6449213 


9994215 


1.6 X 10"" 


5.7 X lO-'' 


C/0.22 


CNTN4' 


rs7652782 


2821616 


3.4 X 10"* 


7.1 X 10-7 


A/0.06 


CNTN4 


rs6786387 


2822150 


2.6 X 10"'* 


7.6 X 10-7 


A/0.03 


CNTN4 


rs6786174 


2982630 


6.1 X 10"'* 


1.2 X 10-'^ 


C/0.08 


CNTN4 


rs 17586876 


2964182 


1.3 X 10-3 


2.0 X 10-'^ 


G/0.09 



"SNf^ Single Nucleotide Polymorphism. 

''MGA, Measured Genotype Analysis. 

'^JLA, Joint linl<age association analysis. 

'^SLC2A9: solute carrier protein 2 family, member 9. 

"CNTNA: Contactin 4. 
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DISCUSSION 

Our GWAS found polymorphisms in SLC2A9 to be signifi- 
cantly associated with SUA levels in Mexican Americans. This 
is the first GWAS for SUA in Mexican Americans and replicates 
results from several studies conducted in European populations 



Table 2 | Association of tlie most significant SNPs in solute carrier 
protein 2 family, member 9 [SLC2A9) gene with serum uric acid 
(mg/dl). 



SNP 


Mean effect 
size (%)= 


lib 


12 


22 


rs6832439 


5.3 


5.6 (1.4)'= 


5.3 (1.4) 


4.8 (1.3) 


rs13131257 


4.9 


5.6 (1.4) 


5.3 (1.4) 


4.8 (1.4) 


rs737267 


5.0 


5.6 (1.4) 


5.3 (1.4) 


4.8 (1.3) 


rs6449213 


4.5 


5.5 (1.4) 


5.2 (1.4) 


4.2 (1.2) 



^Proportion of the residual plienotypic variance tliat is expiained by tlie minor 

allele of the SNR 

^ 1 - major allele; 2 - minor allele. 

'^genotype-specific mean (standard deviation) of serum uric acid levels (mg/dl). 



Q in 
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Genome -wide 






significance level 
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Chromosome 


FIGURE 5 1 A joint linkage-association approacii shows a significant 
LOD score for serum uric acid on chromosome 3. 



from the UK, Germany, Croatia and Sardinia. Recent GWAS 
have shown consistently that single nucleotide polymorphisms 
(SNPs) in the SLC2A9 gene are associated with SUA levels (Li 
et al, 2007; Doring et al, 2008; McArdle et al, 2008; Stark 
et al, 2008; Vitart et al, 2008; Zemunik et al, 2009; Yang 
et al., 2010). A study conducted in a Croatian sample found that 
SLC2A9 variants were associated with SUA with SNPs explain- 
ing 1.7-5.3% of the variance in SUA levels. In the same study, 
this association was replicated in a sample from the island of 
Orkney (Vitart et al., 2008). Two of the three SNPs associated 
with SUA in this study were the same as shown in our study 
(rs6449213 and rs737267). In addition, the effect sizes (propor- 
tion of residual phenotypic variance explained by the SNP) of 
our significant association are in the same range (4-5% of total 
phenotypic variance). Similarly, in a cohort from Germany, SUA 
levels were associated with SLC2A9 variants (19). Of the several 
SNPs that were associated only rs6449213 was common between 
our two studies. They also showed sex-specific effects of vari- 
ants in SLC2A9 on SUA levels. This finding was replicated by 
Brandstatter et al. (2008). Other studies that showed association 
of polymorphisms in SLC2A9 with SUA were conducted in indi- 
viduals from Sardinia (Li et al., 2007), Germany (Stark et al., 

2008) , and the United Kingdom (Wallace et al, 2008; Kolz et al, 

2009) as well as in Asian (Matsuo et al., 2008; Cummings et al., 
2010; Tabara et al, 2010; Tu et al, 2010; Guan et al, 2011) 
and African American populations (Dehghan et al, 2008; Rule 
et al., 2011). To investigate the functional aspects of the SLC2A9 
gene, Caulfield et al. (2008) showed that SLC2A9 was a high- 
capacity urate transporter and can exchange glucose for urate 
in the process of secretion of urate into the urine. They also 
confirmed the previously reported association of SLC2A9 with 
SUA in six different cohorts of European ancestry. Several other 
loci have also been reported to be associated with serum uric 



Table 3 | Solute carrier protein 2 family, member 9 {SLC2A9) SNPs (chromosome 4) associated with cardiovascular or renal disease risk factors. 



Trait 


SNP*" 


MGA" (p-value) 


SNP coordinates 


Minor/major allele 


Minor allele frequency 


Effect size{%)'' 


Body weight 


rs938553 


5.5 X 


10-6 


9925526 


A/G 


0.12 


2.8 


BMI 


rs10003001 


9.8 X 


10-4 


9984475 


A/G 


0.07 


2.1 




rs10008035 


9.1 X 


10-4 


9999335 


A/C 


0.13 


2.0 




rs938553 


7.3 X 


10-4 


9925526 


A/G 


0.12 


2.3 


Waist circunnference 


rs10003001 


5.0 X 


10-4 


9984475 


A/G 


0.07 


2.2 




rs10008035 


8.0 X 


10-4 


9999335 


A/C 


0.13 


2.0 




rs938553 


2.0 X 


10-4 


9925526 


A/G 


0.12 


2.6 


Albumin/creatinine ratio 


rs1014290 


4.0 X 


10-4 


10001861 


G/A 


0.33 


2.1 




rs10805346 


6.0 X 


10-4 


9920347 


A/G 


0.41 


2.0 




rs13129697 


3.0 X 


10-4 


9926967 


C/A 


0.48 


2.2 




rs7660895 


6.0 X 


10-" 


9985445 


A/G 


0.46 


2.0 


Total antioxidant status 


rs6832439 


1.4 X 


10-*^ 


9924319 


A/G 


0.36 


2.9 




rs737267 


1.4 X 


10-6 


9934744 


A/C 


0.38 


2.9 



"Empirical significance based on number of SNPs analyzed in SLC2A9 gene (p < 7.4 x 70"^). 
''SNP Single Nucleotide Polymorphism. 
"MGA, Measured Genotype Analysis. 

'^Proportion of the residual phenotypic variance that is explained by the minor allele of the SNP 
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acid levels in GWAS as well as candidate gene association studies 
(Table 4). 

The solute carrier protein 2 family, member 9 (SLC2A9) gene 
was recently cloned and identified as a member of the solute 
carrier protein 2 family (Le et al., 2008). Originally thought to 
be only a facilitated hexose transporter, it was found to be pri- 
marily involved in uric acid transport. Its two forms SLC2A9a 
and SLC2A9b are expressed in the basolateral and apical mem- 
branes of the kidney while SLC2A9a is also expressed in liver. 
Their amino acid sequences are identical, except that SLC2A9b 
has a shorter and modified N-terminus (Preitner et al., 2009). 
Both forms are involved in renal urate transport (Le et al, 2008; 
Cheeseman, 2009). Caulfield et al. (2008) proposed that extracel- 
lular glucose and intracellular urate are exchanged by the SLC2A9 
transporter. In the kidney SLC2A9 maintains uric acid reab- 
sorption independent of other uric acid transporters (Preitner 
et al., 2009). On one hand SLC2A9 mutations are associated with 
hypouricemia and hyperuricosuria, while on the other, muta- 
tions in SLC2A9 that alter the protein impair urate secretion into 
the urine resulting in hyperuricemia (Le et al., 2008; Cheeseman, 
2009). 

Given the ubiquity of evidence for association of SUA with 
SLC2A9 polymorphisms and the relatively large effect size of 
these associations reported in other studies, it is not surpris- 
ing that we could identify this locus in our study despite its 
modest sample size. In addition, our joint-linkage association 
approach has the potential to maximize the information in a sam- 
ple of related individuals, thus amplifying a signal taking into 



account fixed effects of marker genotypes (association) and the 
random effects of shared sequence identity (linkage). With this 
approach we found suggestive evidence of association on chro- 
mosome 3p25-p26 as we showed in our prior linkage analysis 
(Voruganti et al., 2009a). Associated SNPs were found within the 
contactin 4 (CNTN4) gene, which is located in the one-LOD sup- 
port interval in our previous linkage scan for SUA levels. The 
CNTN4 gene codes for a member of the contactin subgroup of 
cell adhesion molecules of the immunoglobulin (Ig) superfamUy. 
This family has an important role in the formation and func- 
tioning of the neuronal networks; specifically, CNTN4 is known 
to play a key role in development of the central nervous system. 
The possible functional relationship of this gene to SUA levels is 
unknown. 

Our association and JLA results are especially interesting 
given the recent attention to the relative importance of common 
and rare genetic variants for risk of common complex diseases 
(Maher, 2008; Manolio et al., 2009). Common marker variants 
like those in SLC2A9 (minor allele frequency 36-38% in our 
cohort) should, in theory, be in LD with one or more common 
functional variants whose effect on the trait of interest should be 
similar in different lineages. On the other hand, low-frequency 
functional variants (and low- frequency marker alleles in LD with 
them) are expected to appear in only a subset of lineages; these 
pedigree-specific effects are more likely to be identified by linkage 
than by association. In this context, it is interesting that our repli- 
cation of association at SLC2A9 is not matched by a linkage signal, 
while the SNPs identified by joint linkage and association within 



Table 4 | Previously reported association of serum uric acid levels with genes other than SLC2A9. 



SNP 


Joint link- 
age/association 
(JLA) in our 
study (p-value) 


Gene name 


Gene symbol 


Population 


References 


rs2231142 


0.24'' 


ATP-binding cassette family 
G, member 2 


ABCG2 


Mexican Americans, 
African Americans, 
European Americans, 
American Indians 
European Americans 
Japanese 
Whites and Blacks 


Kottgen etal., 2013; 
Zhang etal., 2013 

Kolz etal., 2009 
Tabara et al., 2010 
Dehghan et al., 2008 


rs1471633 


0.79 {rs1967017)'' 


PDZ domain containing 1 


PDZK1 


European Americans 


Kottgen etal., 2013 


rsl2129861 


0.17 {rs1298954)'= 


PDZ domain containing 1 


PDZK1 


European Americans 


Kolz etal., 2009 


rsl171614 


0.11 (rs1171617)'' 


Solute carrier protein family 
16A, member 9 


SLC16A9 


European Americans 


Kolz etal., 2009; 
Kottgen etal., 2013 


rs12800450 


0.05 (rs505802)= 


Solute carrier protein family 
22A, member 12 


SLC22A12 


African Americans 


Tin et al., 2011 


rs2078067 


0.02 

(rs10792438)'= 


Solute carrier protein family 
22A, member 11 


SLC22A11 


European ancestry 


Yang et al., 2010 


rsl 165205 


0.92 (rs9393672)'' 


Solute carrier protein family 
17A, member 3 


SLC17A3 


Whites and Blacks 


Dehghan et al., 2008 



"SNP typed on our array. 
^ Proxy SNP on our array. 

"No proxy SNP available; strongest evidence of JLA in our study annotated to this gene. 

Proxy SNPs identified with SNAP (SNP Annotation and Proxy Search) Version 2.2: https://www.broadlnstltute.org/mpg/snap/ 
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a strong linkage peak on chromosome 3 have low minor allele fre- 
quencies (3-9%) — both findings being consistent with theoretical 
expectation. These results may also explain why our chromosome 
3 QTL has apparently not been reported by other studies, which 
are primarily population-based; however, the current upsurge in 
analysis of next-generation sequencing data and improved meth- 
ods for rare variant detection may provide opportunities to test 
our findings in the near future. 

Since SUA is associated with obesity, type 2 diabetes and CVD 
risk and is advocated by many to be included as a component 
of the metabolic syndrome, we conducted a measured geno- 
type analysis for SLC2A9 SNPs and these risk factors. We found 
evidence of significant associations of SLC2A9 SNPs with body 
weight, BMI and waist circumference and suggestive associations 
with TAS. In the same study we also reported significant genetic 
correlations between SUA and waist circumference in the same 
study (Voruganti et al., 2009a). Variants in the uric acid trans- 
porter gene iSLC2A9) that have been associated with lower SUA 
were also associated with decreased blood pressure (Parsa et al., 
2012). Similarly, another SLC2A9 variant (rsl014290), associ- 
ated with lower SUA were also associated with decreased risk for 
diabetes mellitus (Liu et al., 2011) in Han Chinese. In a study 
by Brandstatter et al. (2008), the association between SLC2A9 
SNPs and SUA was significantly modified by BMI. They showed 
stronger effect size in individuals with higher BMI. However, the 
precise mechanism that may explain this association is not clear. 
In contrast, Vitart et al. (2008) found no association between 
SLC2A9 SNPs and metabolic syndrome components. Caulfield 
et al. (2008) explored the association of SLC2A9 SNPs and blood 
pressure traits and found no significant association. We also 
observed that the SNPs associated with body weight, BMI, or 
waist circumference were different than those associated with 
SUA levels. However, the two SNPs associated with TAS were also 
associated with SUA levels. SUA is known to act as a free rad- 
ical scavenger (Pasalic et al, 2012) and combat oxidative stress 
(lohnson et al., 2011). Therefore, the possible link between them 
may be the role of SUA as an "antioxidant." 

In addition, we found suggestive association between SLC2A9 
SNPs and the renal disease development marker UACR. Knock 
out rodent studies have shown that SLC2A9 knockout mice 
develop urate nephropathy characterized by intra-tubular urate 
lithiasis, interstitial inflammation and fibrosis and tubular atro- 
phy (Preitner et al., 2009). SUA levels in humans are primar- 
ily determined by renal uric acid clearance, with two-thirds of 
uric acid turnover being accounted for by its urinary excre- 
tion. Hyperuricemia is predictably associated with a decrease 
in GFR (CuUeton, 2001; Johnson et al, 2003; Feig, 2009) and 
90% of clinically- recognized hyperuricemia results from impaired 
renal excretion of uric acid (Le et al., 2008). Abnormal renal 
urate handling can promote both uric acid and calcium oxalate 
nephrolithiasis by distinct mechanisms; low urine pH favors uric 
acid crystallization while mixed uric acid- or sodium urate- 
calcium oxalate stones develop by epitaxy, "salting out" of cal- 
cium oxalate or adsorption of crystal growth inhibitors (Liebman 
et al., 2007). Hyperuricemia has also been implicated in the 
development of hypertension, a major risk factor for the devel- 
opment and progression of chronic kidney disease (Feig, 2009). 



Micropuncture studies have shown that elevated uric acid levels 
cause cortical vasoconstriction and increase glomerular capil- 
lary pressure leading to progressive glomerular injury (Sanchez- 
Lozada et al., 2002). Given the role of SLC2A9 in renal urate 
transport, the association between SLC2A9 SNPs and UACR 
assumes pathogenic significance. 

To summarize, this is the first study to report the association of 
polymorphisms in the SLC2A9 gene with SUA levels in a genome- 
wide analyses in Mexican Americans. In addition, this is first 
time that any association between SLC2A9 SNPs and BMI, body 
weight, waist circumference, UACR, and TAS has been reported. 

AUTHOR CONTRIBUTIONS 

Conceived and designed the experiments: Venkata Saroja 
Voruganti, Jack W. Kent Jr., John Blangero, Anthony G. 
Comuzzie; Performed or supervised all aspects of statistical anal- 
yses: Venkata Saroja Voruganti, Jack W. Kent Jr., John Blangero, 
Anthony G. Comuzzie; Helped with statistical analyses: Harald H. 
H. Goring, Laura Almasy, Thomas D. Dyer, Michael C. Mahaney, 
John Blangero, Anthony G. Comuzzie; Helped with prepara- 
tion/editing of manuscript: Subrata Debnath, Shelley A. Cole, 
Karin Haack, Harald H. H.Goring, Melanie A. Carless, Joanne 
E. Curran, Matthew P. Johnson, Laura Almasy, Thomas D. Dyer, 
Jean W. MacCluer, Eric K. Moses, Hanna E. Abboud, Michael 
C. Mahaney, John Blangero, Anthony G. Comuzzie; Wrote the 
manuscript: Venkata Saroja Voruganti, Jack W. Kent Jr. 

ACKNOWLEDGMENTS 

We wish to thank all participants of SAFHS for their cooperation 
and generous participation. Our work was supported by grant 
POl HL045522 from National Heart, Lung and Blood Institute 
(NHLBI) and grant MH059490 from National Institute of Mental 
Health (NIMH) of the National Institutes of Health. A dona- 
tion from the Azar and Shepperd families in San Antonio paid 
for the SNP typing in the SAFHS. This investigation was con- 
ducted in part in facilities constructed with support from the 
Research Facilities Improvement Program under grant numbers 
C06 RR013556 and C06 RR017515. The AT&T Genomics com- 
puting center supercomputing facilities used for this work were 
supported in part by a gift from the AT&T Foundation and with 
support from the National Center for Research Resources Grant 
Number S10RR029392. The funders had no role in study design, 
data collection and analysis, decision to publish, or preparation 
of the manuscript. 

REFERENCES 

Abecasis, G. R., Cherny, S. S., Cookson, W. O., and Cardon, L. R. (2002). Merlin- 
rapid analysis of dense genetic maps using sparse gene flow trees. Nat. Genet. 30, 
97-101. doi: 10.1038/ng786 

Almasy, L., and Blangero, J. (1998). Multipoint quantitative-trait linkage analysis in 
general pedigrees. Am. }. Hum. Genet 62, 1198-1211. doi: 10.1086/301844 

Arar, N. H., Voruganti, V. S., Nath, S. D., Thameem, E, Bauer, R., Cole, S. 
A., et al. (2008). A genome-wide search for linkage to chronic kidney dis- 
ease in a community-based sample: the SAFHS. Nephrol Dial. Transplant 23, 
3184-3191. doi: 10.1093/ndt/gfn215 

Boerwinkle, E., Chakroborty, R., and Sing, C. F. (1986). The use of measured geno- 
type information in the analysis of quantitative phenotypes in man. I. Models 
and analytical methods. Ann. Hum. Genet. 50, 181-194. doi: 10.1111/j.l469- 
1809.1986.tb01037.x 



www.frontiersin.org 



December 2013 | Volume 4 | Article 279 | 7 



Voruganti et al. 



SUA GWAS in Mexican Americans 



Brandstatter, A., Kiechl, S., Kollerits, B., Hunt, S. C, Heid, I. M., Coassin, S., et al. 

(2008). Sex-specific association of the putative fructose transporter SLC2A9 

variants with uric acid levels is modified by BMI. Diabetes Care 31, 1662-1667. 

doi: 10.2337/dc08-0349 
Cameron, J. S., and Simmonds, H. A. (2005). Hereditary hyperuricemia and renal 

disease. Semin. Nephrol. 25, 9-18. doi: 10.1016/j.semnephrol.2004.09.003 
Caulfield, M. )., Munroe, P. B., O'Neill, D., Wilkowska, K., Charchar, F. J., Doblado, 

M., et al. (2008). SLC2A9 is a high-capacity urate transporters in humans. PLoS 

Med. 5:el97. doi: 10.1371/journal.pmed.0050197 
Charles, B. A., Shriner, D., Doumatey, A., Chen, G., Zhou, )., Huang, H., et al. 

(2011). A genome-wide association study of serum uric acid in African 

Americans. £MD Med. Genomics 4, 4-17 . doi: 10.1186/1755-8794-4-17 
Cheeseman, C. (2009). Solute carrier family 2, member 9 and uric 

acid homeostasis. Curr. Opin. Nephrol. Hypertens. 18, 428^32. doi: 

10.1097/MNH.0b013e32832ee3de 
Cirillo, P., Sato, W., Reungjui, S., Heinig, M., Gesch, M., Saulin, Y., et al. (2006). 

Uric acid, the metabolic syndrome, and renal disease. /. Am. Soc. Nephrol. 17, 

S165-S168. doi: 10.1681/ASN.2006080909 
Culleton, B. F. (2001). Uric acid and cardiovascular disease: a renal-cardiac rela- 
tionship? Cwrr. Opin. Nephrol. Hypertens. 10, 371-375. doi: 10.1097/00041552- 

200105000-00012 

Citmmings, N., Dyer, T. D., Kotea, N., Kowlessur, S., Chitson, P., Zimmet, P., et al. 
(2010). Genome-wide scan identifies a quantitative trait locus at 4pl5.3 for 
serum urate. Eur. J. Hum. Genet. 18, 1243-1247. doi: 10.1038/ejhg.2010.97 

Dehghan, A., Kottgen, A., Yang, Q., Hwang, S. J., Kao, W. L., Rivadeneira, F., 
et al. (2008). Association of three genetic loci with uric acid concentration 
and risk of gout: a genome-wide association study. Lancet 372, 1953—1961. doi: 
10.1016/50140-6736(08)61343-4 

Dixon, A. S. (1960). Famihal hyperuricemic nephropathy. Proc. R. Soc. Med. 53, 
967-968. 

Domagk, G. R, and Schlicke, H. H. (1968). A colorimetric method using uricase 

and peroxidase for the determination of uric acid. Anal. Biochem. 22, 219-224. 

doi: 10.1016/0003-2697(68)90309-6 
Doria, A., and Krolewski, A. S. (2011). Lowering' serum uric acid levels to prevent 

kidney failure. Nat. Rev. Nephrol. 7, 495-496. 
Dbring, A., Gieger, C, Mehta, D, Gohlke, H., Prokisch, H., Coassin, S., et al. (2008). 

SLC2A9 influences uric acid concentrations with pronounced sex-specific 

effects. Nat Genet 40, 430^36. doi: 10.1038/ng.l07 
Ejaz, A. A., Dass, B., Kambhampati, G., Ejaz, N. I., Maroz, N., Dhatt, G. S., 

et al. (2012). Lowering serum uric acid to prevent acute kidney injury. Med. 

Hypotheses 78, 796-799. doi: 10.1016/i.mehy2012.03.011 
Feig, D. I. (2009). Uric acid: a novel mediator and marker of risk in 

chronic kidney disease? Curr. Opin. Nephrol. Hypertens. 18, 526-530. doi: 

10.1097/MNH.Ob013e328330d9dO 
Flegal, K. M., Ezzati, T. M., Harris, M. L, Haynes, S. G., Juarez, R. Z., Knowler, 

W. C, et al. (1991). Prevalence of diabetes in Mexican Americans, Cubans, and 

Puerto Ricans from the Hispanic Health and Nutrition Examination Survey, 

1982-1984. Diabetes Care.U, 628-638. doi: 10.2337/diacare.l4.7.628 
Friedlander, Y., Kark, J. D., and Stein, Y. (1988). Family resemblance for serum 

uric acid in a Jerusalem sample of families. Hum. Genet. 79, 58-63. doi: 

10.1007/BF00291711 

Go, A. S. (2013). Heart disease and stroke statistics-2013 update: a 
report from the American Heart Association. Circulation 127, e6. doi: 
10. 1 161/CIR.ObO 1 3e3 1 828 1 24ad 

Gotsman, I., Kereb, A., Lotan, C, and Zwas, D. R. (2012). Changes in uric acid 
levels and allopurinol use' in chronic heart failure: association with improved 
survival./. Card. Fail. 18, 694-701. doi: 10.1016/j.cardfail.2012.06.528 

Guan, M., Zhou, D., Ma, W., Chen, Y., Zhang, J., and Zou, H. (2011). Association 
of an intronic SNP of SLC2A9 gene with serum uric acid levels in the Chinese 
male Han population by high-resolution melting method. Clin. Rheumatol. 30, 
29-35. doi: 10.1007/sl0067-010-1597-x 

Johnson, R. J., Kang, D. H., Feig, D., Kivlighn, S., Kanellis, J., Watanabe, 
S., et al. (2003). Is there a pathogenetic role for uric acid in hyperten- 
sion and cardiovascular and renal disease? Hypertension 41, 1183-1190. doi: 
10.1161/01.LiYR0000069700.62727.C5 

Johnson, R. J., Lanaspa, M. A., and Gaucher, E. A. (2011). Uric acid: a danger 
signal from the RNA world that may have a rt)le in the epidemic of obe- 
sity, metabolic syndrome, and cardiorenal disease: evolutionary considerations. 
Semin. Nephrol. 31, 394-399. doi: 10.1016/i.semnephrol.2011.08.002 



Karns, R., Zhang, G., Sun, G., Rao Indugula, S., Cheng, H., Havas-Augustin, D., 
et al. (2012). Genome-wide association of serum uric acid concentration: reph- 
cation of sequence variants in an island populations of the Adriatic coast of 
Croatia. Ann. Hum. Genet. 75, 121-127. doi: 10.1111/j.l469-1809.2011.00698.x 

Kolz, M., Johnson, T., Sanna, S., Teumer, A., Vitart, V., Perola, M., et al. (2009). 
Meta-analysis of 28,141 individuals identifies common variants within five 
new loci that influence uric acid concentrations. PLoS Genet. 5:el000504. doi: 
10.1371/journal.pgen.l000504 

Kottgen, A., Albrecht, E., Teumer, A., Vitart, V., Krumsiek, J., Hundertmark, C, 
et al. (2013). Genome-wide association analyses identify 18 new loci associ- 
ated with serum urate concentrations. Nat. Genet. 45, 145-154. doi: 10.1038/ 
ng.2500 

Le, M. T., Shafiu, M., Mu, W., and lohnson, R. J. (2008). SLC2A9-a fructose trans- 
porter identified as a novel uric acid transporter. Nephrol. Dial. Transplant. 23, 
2746-2749. doi: 10.1093/ndt/gfn349 

Li, S., Sanna, S., Maschio, A., Busonero, F., Usala, G., Mulas, A., et al. (2007). The 
GLUT9 gene is associated with serum uric acid levels in Sardinia and Chianti 
cohorts. PLoS Genet 3:el94. doi: 10.1371/journal.pgen.0030194 

Liebman, S. E., Taylor, J. G., and Bushinsky, D. A. (2007). Uric acid nephrolithiasis. 
Curr. Rheumatol. Rep. 9, 251-257. doi: 10.1007/sll926-007-0040-z 

Liu, W. C, Hung, C. C, Chen, S. C, Lin, M. Y, Chen, L. I., Hwang, D. Y, et al. 
(2011). The rsl014290 polymorphisms of the SLC2A9 gene is associated with 
type 2 diabetes mellitus in Han Chinese. Exp. Diabetes Res. 2011:527520. doi: 
10.1155/2011/527520 

MacCluer, J. W., Stern, M. P., Almasy, L., Atwood, L. A., Blangero, J., Comuzzie, A. 
G., et al. (1999). Genetics of atherosclerosis risk factors in Mexican Americans. 
Nutr.Rev. 57, S59-S65. doi: 10.1111/j.l753-4887.1999.tb01790.x 

Maher, B. (2008). Personal genomes: the case of the missing heritability. Nature 
456, 18-21. doi: 10.1038/456018a 

Manolio, T. A., Collins, E S., Cox, N. J., Goldstein, D. B., Hindorff, L. A., Hunter, 
D. J., et al. (2009). Finding the missing heritability of complex diseases. Nature 
461, 747-753. doi: 10.1038/nature08494 

Matsuo, H., Chiba, T., Nagamori, S., Nakayama, A., Domoto, H., Phetdee, K., 
et al. (2008). Mutations in glucose transporter 9 gene SLC2A9 cause renal 
hypouricemia. Am. J. Hum. Genet. 83, 744-751. doi: 10.1016/j.ajhg.2008. 
11.001 

McArdle, R E, Parsa, A., Chang, Y. R, Weir, M. R., O'Connell, J. R., Mitchell, B. D, 
et al. (2008). Association of a common nonsynonymous variant in GLUT9 with 
serum uric acid levels in old order amish. Arthritis Rheum. 58, 2874-2881. doi: 
10.1002/art.23752 

Mitchell, B. D., Almasy, L. A., Rainwater, D. L., Schneider, J. L., Blangero, J., Stern, 
M. P., et al. (1999). Diabetes and hypertension in Mexican American fam- 
ilies: relation to cardiovascular risk. Am. J. Epidemiol. 149, 1047-1056. doi: 
10.1093/oxfordjournals.aje.a009750 

Nakagawa, T., Kang, D.-H., Feig, D., Sanchez-Lozada, L. G., Srinivas, T. R., 
Ejaz, A. A., et al. (2006). Unearthing uric acid: an ancient factor with 
recently found significance in renal and CVD. Kidney Int. 69, 1722-1725. doi: 
10.1038/sj.ki.5000391 

Nakagawa, T., Tuttle, K. R., Short, R. A., and Johnson, R. J. (2005). Hypothesis: 
fructoseinduced hyperuricemia as a causal mechanism for the epidemic of the 
metabolic syndrome. Nat. Clin. Pract. Nephrol. 1, 80-86. doi: 10.1038/ncp- 
neph0019 

Nakayama, A., Matsuo, H., Takada, T., Ichida, K., Nakamura, T., Ikebuchi, Y, et al. 

(2011). ABCG2 is a high-capacity urate transporter and its genetic impairment 

increases serum uric acid levels in humans. Nucleosides Nucleotides Nucleic Acids 

30, 1091-1097. doi: 10.1080/15257770.2011.633953 
Nath, S. D., Voruganti, V. S., Arar, N. H., Thameem, F., Lopez-Alvarenga, J. C, 

Bauer, R., et al. (2007). Genome scan for determinants of serum uric acid 

variability./. Am. Soc. Nephrol. 18, 3156-3163. doi: 10.1681/ASN.2007040426 
Parsa, A., Brown, E., Weir, M. R., Fink, J. C, Shuldiner, A. R., Mitchell, B. D., 

et al. (2012). Genotype-based changes in serum uric acid affect blood pressure. 

Kidney Int 81, 502-507. doi: 10.1038/ki.2011.414 
Pasalic, D., Marinkovic, N., and Feher-Turkovic, L. (2012). Uric acid as one of the 

important factors in multifactorial disorders: facts and controversies. Biochem. 

Med. (Zagrab) 22, 63-75. doi: 10.11613/BM.2012.007 
Preitner, E, Bonny, O., Laverriere, A., Rotman, S., Firsov, D., Da Costa, A., et al. 

(2009). Glut9 is a major regulator of urate homeostasis and its genetic inac- 

tivation induces hyperuricosuria and urate nephropathy. Proc. Natl. Acad. Sci. 

U.S.A. 106, 15501-15506. doi: 10.1073/pnas.0904411106 



Frontiers in Genetics | Applied Genetic Epidenniology 



December 2013 | Volume 4 | Article 279 | 8 



Voruganti et al. 



SUA GWAS in Mexican Americans 



Rajendra, N. S., Ireland, S., George, J., Beich, J. J., Lang, C. C, and Struthers, A. D. 
(2011). Mechanistic insights into the therapeutic use of high-dose allopurinol 
manqlna pectoris. /. Am. Coli. Cardiol. 58, 820-828. doi: 10.1016/j.jacc.2010. 
12.052 

Rao, D. C, Laskarzewski, R M., Morrison, J. A., Kelly, K., and Glueck, C. J. (1982). 
The clinical Lipid Research Clinic Family Study: familial determinants of plasma 
uric acid. Hum. Genet. 60, 257-261. doi: 10.1007/BF00303013 

Rice, T., Vogler, G. R, Perry, T. S., Laskarzewski, R M., Province, M. A., and Rao, D. 
C. (1990). Heterogenity in the familial aggregation of fasting uric acid levels in 
five North American populations: the Lipid Research Clinics Family Study. Am. 
J. Med. Genet. 36, 219-225. doi: 10.1002/ajmg.l320360216 

Rule, A. D., de Andrade, M., Matsumoto, M., Mosley, T. H., Kardia, S., and 
Turner, S. T. (2011). Association between SLC2A9 transporter gene variants and 
uric acid phenotypes in African. American and white families. Rheumatology 
(Oxford) 50, 871-878. doi: 10.1093/rheumatology/keq425 

Rule, A. D., Fridley, B. L., Hunt, S. C, Asmann, Y., Boerwinkle, E., Pankow, 
J. S., et al. (2009). Genome-wide linkage analysis for uric acid in fami- 
Ues enriched for hypertension. Nephrol. Dial. Transplant. 24, 2414-2420. doi: 
10.1093/ndt/gfp080 

Sanchez-Lozada, L. G., Tapia, E., Avila-Casado, C, Soto, V., Franco, M., Santamaria, 
J., et al. (2002). MUd hyperuricemia induces glomerular hypertension in nor- 
mal rats. Am. }. Physiol. Renal Physiol. 283, F1105-F1110. doi: 10.1152/ajpre- 
nal.00170.2002 

Self, S. G., and Liang, K.-Y. (1987). Asymptotic properties of maximum likelihood 

estimators and likelihood ratio tests under nonstandard conditions. /. Am. Stat. 

Assoc. 82, 605-610. doi: 10.1080/01621459.1987.10478472 
Sobel, E., and Lange, K. (1996). Descent graphs in pedigree analysis: Applications 

to haplotyping, location scores, and marker-sharing statistics. Ann. Hum. Genet. 

58, 1323-1337. 

Stark, K., Reinhard, W., Neureuther, K., Wiedmann, S., Sedlacek, K., Baessler, A., 
et al. (2008). Association of common polymorphisms in GLUT9 gene with gout 
but not with coronary artery disease in a large case-control study. PLoS ONE 
3:el948. doi: 10.1371/journaLpone.0001948 

Tabara, Y., Kohara, K., Kawamoto, R., Hiura, Y., Nishimura, K., Morisaki, T, 
et al. (2010). Association of four genetic loci with uric acid levels and reduced 
renal function: the J-SHIPP Suita study Am. J. Nephrol. 32, 279-286. doi: 
10.1159/000318943 

Tang, W., Miller, M. B., Rich, S. S., North, K. E., Pankow, J. S., Borecki, I. B., 
et al. (2003). Linkage analysis of a composite factor for the multiple metabolic 
syndrome. The National Heart, Lung and Blood Institute Family Heart Study. 
Diabetes 52, 2840-2847. doi: 10.2337/diabetes.52.11.2840 

Tin, A., Woodward, O. M., Kao, W. H., Liu, C. X, Lu, X., Nails, M. A., et al. 
(2011). Genome-wide association study for serum urate concentrations and 
gout among African Americans identifies genomic risk loci and a novel URAT 
1 loss-of-fUnction allele. Hum. Mot Genet 20, 4056-4068. doi: 10.1093/ 
hmg/ddr307 

Tu, H. P, Chen, C. J., Tovosia, S., Ko, A. M., Lee, C. H., Ou, T. T, et al. (2010). 
Associations of a non-synonymous variant in SLC2A9 with gouty arthritis and 
uric acid levels in Han Chinese subjects and Solomon Islanders. Ann. Rheum. 
Bis. 69, 887-890. doi: 10.1 136/ard.2009.1 13357 

Vitart, v., Rudan, I., Hayward, C, Gray N. K., Floyd, )., Palmer, C. N., et al. 
(2008). SLC2A9 is a newly identified urate transporter influencing serum 
urate concentration, urate excretion and gout. Nat. Genet. 40, 437-442. doi: 
10.1038/ng.l06 

Voruganti, V. S., Nath, S. D., Cole, S. A., Thameem, F, Jowett, J. B., Bauer, R., 
et al. (2009a). Genetics of variation in serum uric acid and cardiovascular risk 



factors in Mexican-Americans. /. Clin. Endocrinol. Metab. 94, 632-638. doi: 
10.1210/jc.2008-0682 

Voruganti, V. S., Goring, H. H. H., Mottl, A., Franceschini, N., Haack, K., Laston, 
S., et al. (2009b). Genetic influence on variation in serirm uric acid in American 
Indians: the Strong Heart Family Study. Hum. Genet. 126, 667-676. doi: 
10.1007/s00439-009-0716-8 

Wallace, C, Newhouse, S. ]., Braund, P., Zhang, E, Tobin, M., Ealchi, M., et al. 
(2008). Genome-wide association study identifies genes for biomarkers of car- 
diovascular disease: serum urate and dyslipidemia. Am. ]. Hum. Genet. 82, 
139-149. doi: 10.1016/j.ajhg.2007.11.001 

Wilk, J. B., Djousse, L., Borecki, I., Atwood, L. D., Hunt, S. C, Rich, S. S., et al. 
(2000). Segregation analysis of serum uric acid in the NHLBI Family Heart 
Study Hum. Genet. 106, 355-359. doi: 10.1007/s004390051050 

Woodward, O. M., Kottgen, A., Coresh, J., Boerwinkle, E., Guggino, W. B., and 
Kottgen, M. (2009). Identification of a urate transporter, ABCG2, with a com- 
mon functional polymorphism causing gout. Proc. Natl. Acad. Sci. U.S.A. 106, 
10338-10342. doi: 10.1073/pnas.0901249106 

Yang, Q., Guo, C. Y., Cupples, A., Levy, D., Wilson, R W., and Fox, C. S. (2005). 
Genome-wide search for genes affecting serum uric acid levels: the Framingham 
Heart Study Metabolism 54, 1435-1441. doi: 10.1016/j.metabol.2005.05.007 

Yang, Q., Kottgen, A., Dehghan, A., Smith, A. V., Glazer, N. L., Chen, M. H., et al. 
(2010). Multiple genetic loci influence serum urate levels and their relation- 
ship with gout and cardiovascidar disease risk factors. Circ. Cardiovasc. Genet. 
3, 523-530. doi: 10.1161/CIRCGENETICS.109.934455 

Yelken, B., Caliskan, Y, Gorgulu, N., Altun, I., Yimaz, A., Yazici, H., et al. (2012). 
Reduction of uric acid levels with allopurinol treatment improves endothelial 
function in patients with chronic kidney disease. Clin. Nephrol. 11, 275—282. 
doi: 10.5414/CN107352 

Zemunik, T., Boban, M., Lauc, G., Jankovic, S., Rotim, K., Vatavuk, Z., et al. (2009). 
Genome-wide association study of biochemical traits in Korcula Island, Croatia. 
Croat. Med. J. 50, 23-33. doi: 10.3325/cmj.2009.50.23 

Zhang, L., Spencer, K. L., Voruganti, V. S., Jorgensen, N. W., Fornage, M., Best, L. 
G., et al. (2013). Association of functional polymorphisms rs2231I42 (Q141K) 
in the ABCG2 gene with serum uric acid and gout in 4 US populations: the 
PAGE study Am. J. Epidemiol, doi: 10.1093/aje/kws330. [Epub ahead of print]. 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 26 August 2013; paper pending published: 24 September 2013; accepted: 23 
November 2013; published online: 16 December 2013. 

Citation: Voruganti VS, Kent JW Jr., Debnath S, Cole SA, Haack K, Goring HHH, 
Carless MA, Curran JE, Johnson MP, Almasy L, Dyer TO, MacCluer JW, Moses EK, 
Abboud HE, Mahaney MC, Blangero J and Comuzzie AG (2013) Genome-wide asso- 
ciation analysis confirms and extends the association of SLC2A9 with serum uric acid 
levels to Mexican Americans. Front. Genet. 4:279. doi: 10.3389/fgene.2013.00279 
This article was submitted to Applied Genetic Epidemiology, a section of the journal 
Frontiers in Genetics. 

Copyright © 2013 Voruganti, Kent, Debnath, Cole, Haack, Goring, Carless, Curran, 
Johnson, Almasy, Dyer, MacCluer, Moses, Abboud, Mahaney, Blangero and Comuzzie. 
This is an open-access article distributed under the terms of the Creative Commons 
Attribution License (CC BY). The use, distribution or reproduction in other forums is 
permitted, provided the original author(s) or licensor are credited and that the original 
publication in this journal is cited, in accordance with accepted academic practice. No 
use, distribution or reproduction is permitted which does not comply with these terms. 



www.frontiersin.org 



December 2013 | Volume 4 | Article 279 | 9 



